00010	INTERNAL PROCEDURE LPC;
00020	BEGIN "LPC"
00030	DEFINE ⊂="COMMENT"; ⊂ 10/9/73;
00040	DEFINE ⊃=" "; ⊂ DEFINE ⊃="SAFE"; ⊂ ALTERNATE DEFINITION WHEN SAFE;
00050	REQUIRE "LPC1[X,ALS]" LOAD_MODULE;
00060	FORTRAN REAL PROCEDURE SQRT(REAL X);
00070	FORTRAN REAL PROCEDURE ALOG10(REAL X);
00080	FORTRAN REAL PROCEDURE COS(REAL X);
00090	FORTRAN REAL PROCEDURE SIN(REAL X);
00100	INTEGER DPPOINT,DPP1,DPP2,DATSHIFT;
00110	EXTERNAL FORTRAN PROCEDURE LPC1
00120	    (REFERENCE REAL A,B,R0,C;REFERENCE INTEGER N,I,J);
00130	DEFINE BYTE="((ILDB(BPT) LSH 24)%2↑24)";
00140	⊃ INTERNAL INTEGER ARRAY LIST[0:1];	⊂ Needed but not used;
00150	⊃ INTERNAL REAL ARRAY A,B,C[0:256];
00160	REAL X,SX; ⊃ REAL ARRAY WINDOW[0:256];
00170	INTERNAL REAL R0 ;
00180	INTEGER I,J,L,M,N,BPT,BYTE;
00190	LABEL START;
00200	 M←8; N←2↑M; NF←2*N; DEFINE PI="3.141592653";
00210	FOR I←0 STEP 1 UNTIL N DO  WINDOW[I]←(1-COS((2*PI*I)/N))/2;
00220	FOR I←0 STEP 1 UNTIL N-1 DO A[I]←BYTE;
00230	 FOR I←0 STEP 1 UNTIL N-2 DO A[I]←(A[I+1]-A[I])*WINDOW[I];
00240	 ⊂  LOADS DATA IN A, DIFFERENTIATES AND WINDOWS ;
00250	I←24; J←N%2; LPC1(A[0],B[0],R0,C[0],N,I,J);